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ABSTRACT 

We report results of a theoretical investigation of polarization of the X-ray emissions induced in charge- 
exchange collisions of fully stripped solar wind ions C^"^ and O^"^ with the heliospheric hydrogen atoms. The 
polarization of X-ray emissions has been computed for line-of-sight observations within the ecliptic plane 
as a function of solar wind ion velocities, including a range of velocities corresponding to the slow and fast 
solar wind, and Coronal Mass Ejections. To determine the variability of polarization of heliospheric X-ray 
emissions, the polarization has been computed for solar minimum conditions with self-consistent parameters 
of the solar wind plasma and heliospheric gas and compared with the polarization calculated for an averaged 
solar activity. We predict the polarization of charge-exchange X-rays to be between 3% and 8%, depending on 
the line-of-sight geometry, solar wind ion velocity, and the selected emission lines. 

Subject headings: atomic processes - interplanetary medium - X-rays: diffuse background - X-rays: polariza- 
tion 



1. EVITRODUCTION 

Charge-exchange (CX) collisions between highly charged 
Solar Wind (SW) ions and neutral gas were recently identi- 
fied as an efficient mechanism for production of EUV a nd soft 
X-ray emissions (Robertson et al. 2001, Koutroumpa et al. 



2006] |2009b[ [Bhardwaj et al. 2007, .Snowden„2009) . There 
are indications that the CX collisions in the Heliospliere and 
Geocorona yield between 50% and 80% of the observed soft 
(below 1 keV) X-ray photons, making them a significant con- 
tributor to the soft X-ray background. Heliospheric CX X-ray 
emissions are sensitive to the parameters of the SW plasma 
and strong correlations between variations in the SW inten- 
sity and composition, and intensity of the soft X-ray back- 
ground were observed and anal yzed ( [Robertson et al.|[200H 
|2008[|Koutroumpa et al.|2 009b ). These findings indicate that 
the heUospheric CX X-rays could be used for diagnostics of 
the solar wind composition and velocities, as well as an inde- 
pendent probe of spatial distribution of the heliospheric neu- 
tral gas. The X-ray emission, if measured, could provide ad- 
ditional insight into interaction of the SW plasma with the 
neutral heliospheric gas. Parameters of this interaction are 
critical for modeling of the interstellar gas distribution and 
for interpretation of observational dat a on the diffuse soft X- 
ray emission from the Local Bubble (Lallement et al.]|2005 



[Henley & Shelton|2008l [Welsh & Shelton|jOO^K 

While polarization or optical emissions is readily used in 
investigation of various astrophysical objects, detection of X- 
ray polarization remains a technically challenging feat. Cur- 
rently available observational data do not include polariza- 
tion. However, the angular di stribution of polarized emissions 
is known to be anisotropic ( [Fano & Macek[[l973t [Tanuma[ 
[et al.[[2000| l. If the CX X-ray emissions are strongly polar- 
ized, the total intensities measured in current satellite obser- 
vations should be redefined to adjust for the anisotropy. It 
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is realistic to expect that X-ray polarization will be investi- 
gated in future space missions, making an accurate theoret- 
ical consideration of these phenomena necessary. Currently, 
polarization data are available from laboratory experiments 
for some CX collisions of atoms and ions of astrophysical in- 
terest. For example, polarization spectroscopy of 0^^(Ls'^3/:') 
produced in collisions of O^"*" with He and H2, for projectile 
velocities from 740 to 1200 km s"', showed that the CX X- 
rays are p olarized and stron gly dependent on the projectile 
velocity (Tanuma et al. 200()). We expect CX X-ray emis- 
sions produced in astrophysical environments, su ch as Jupiter 
heavy ion auroras ( Kharchenko et al. 2006 2008| l and coronal 
mass ejections, to exhibit similar properties. 

In this Letter, we report the results of our theoretical study 
of polarization of the X-rays emitted in charge-exchange col- 
lisions between fully stripped SW ions and neutral helio- 
spheric hydrogen. In particular, we investigate the depen- 
dence of the X-ray polarization on distribution of the SW 
plasma, heliospheric neutral gas and SW ion velocity. 

2. POLARIZATION OF THE EMISSIONS INDUCED IN 
CHARGE-EXCHANGE COLLISIONS 

Descriptions of a theoretical model of charge-exchange- 
induced X-rays in the Solar System and calculat ed emission 
spectra have been previously published (Kharchenk o et al. 
1998 2008 ; Kharchenko & Dalgarno 2000 ; Koutroump TetaE 
|2006, |Bhardwaj et al.|2007 ). We consider collisions between 
heavy SW ions and hydrogen gas, X^^ + H ^ X*(2-i)+ + H+, 
taking into account that H is a predominant component of in- 
terstellar gas. This model can be easily extended to include 
minor components of the interstellar gas such as He. 

The line-of-sight (LOS) intensity, expressed as number of 
photons per square centimeter per second, of an emission line 
of the wavelength A is given by 
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where Nxo+ir) is the density of solar wind ions, nu(r) is the 
density of hydrogen, Vrei is the relative collision velocity, 
(Tjj xe* '^he electron capture cross section for an excitation 
of the X-ray emission with the wavelength A. In Eq. ([T]l, Rd 
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is the position of the X-ray detector, and R\i is the distance to 
the heliopause, as defined for a particular geometry and distri- 
bution of heliospheric plasma. For simplicity, we assume that 
Vjei w V, where v is the velocity of SW ions. The polarization 
of the heliospheric X-ray emission in LOS observations can 
be defined as P = (/|| where /|| and /j^ are the 

detected intensities, defined below. 

Projections of angular momenta of the electronic states of 
the excited SW ions produced in CX collisions are oriented 
along the local velocity direction of the SW plasma. Assum- 
ing that a satellite detector can analyze linear polarization of 
X-ray emissions, the intensity /(■(/;) of the radiation, emitted 
from a selected small volume of the heliosphere and normal- 
ized to a sin gle ion emission, can be expressed as (Fano & 
|Macek| 1973) 1 

^(^) = (l - {h'^'Ar^lh^'^AfM)) , (2) 

where R is the distance between the detector and the CX event, 
and h^^^ is the ratio of recoupling coefficients which contain 
the elements of orthogon al transformations betw een the initial 
and final quantum state ( Fano & Macek|1973| l. The angle tp 
indicates the orientation of the linear polarization analyzer. 
In the assumed geometry, we define the polarization parallel 
and perpendicular to the ecliptic plane as /|| = I{tp = 0) and 
I± = = TT /2). Here, AjJ*^* and Aj^' are alignment parameters 
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A^^'(V;)= -A§°' sin^ 61 COS 2V', 



(3) 



is the alignment tensor, 9 is the angle between the 
ion propagation direction and the detection axis. For the de- 
fined geometry, the alignment tensor Aq°' is a scalar which can 
be expressed in terms of the partial cross section a{mi,v) for 
electron capture into the atomic state |y',m,): 



A!^°' = 



(4) 



The parallel and perpendicular intensities in the detector 
frame are then 



I± = 



(5) 
(6) 



Since the cross sec tions a(mi,v) depe nd on the relative veloc- 
ity of collisions v (| Harel et al.|1998| l, the polarization is also 
velocity dependent. 

In order to calculate the X-ray polarization along a LOS, 
we adopt planar geometry with Earth (the X-ray detector is 
assumed to be on a satellite in Earth's orbit). Sun and the CX 
event located in the ecliptic plane and inside the heliosphere. 
We consider Earth to be approximately at the Fall point (Rq = 
1 AU) and specify the LOS direction by the angle a, such 
that for a = the LOS points toward Sun, for a = 7r/2 in the 
direction opposite to the interstellar flow (upwind), and for 
a = -7r/2 along the interstellar flow (downwind). Here, R and 
r are the distances from Earth (detector) and Sun to the CX 
event, respectively. 



We first illustrate general features and estimate the X-ray 
polarization using simplified heliospheric distributions of the 
SW plasma and neutral gas. If the heliospheric hydrogen den- 
sity is assumed to be constant, nuir) = n[^'\ and distribution 
of SW ions to be isotropic as seen from the Sun's location, 

NxQ+(r) = N^^+ (Rq /r) ^, where N^g^ is the ion density at 1 AU, 
we can analytically integrate Eqs. (|6]l up to the heliopause, 
which we assume to be a sphere of the radius R^. 
Resulting expressions for the LOS angle < a < tt are 
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where D = A^^^j+neVfTH.xe+^o/STr, and /? = 

arcsin [(/?()//?/,) sin a] . The X-ray intensities are given 

in cm"^ s"'. The simplified analytic model was used for 
comparison with numerical solutions obtained in the next 
section. It overestimates the X-ray intensities and degree of 
polarization by about a factor three when compared to more 
realistic gas and plasma distributions used in the next section. 

3. POLARIZATION FOR REALISTIC DISTRIBUTIONS OF THE 
HELIOSPHERIC GAS AND PLASMA 

We calculated polarization of X-rays for single-photon re- 
laxation from the 4p, 5p and 6p excited states of C**^^^^ and 
Q*(7+) produced in CX collisions of C^"^ and O^"^ ions wifli 
heliospheric hydrogen. Probability to populate these states 
is an order of magnitude or more h igher than the pr obabil- 
ity to populate other excited states ( Harel et al.| 1998), and 
the selected ions characterize the slow SW, wh ich dominates 
the ecliptic plane ( Koutroumpa et al.|2006 |2009a). Radiative 
cascade from the excited states was not incluHed in the calcu- 
lation since it does not contribute significantly to populations 
of considered excited states. Electron capture cross sections 
depend on the relative velocity of colliding particles; we com- 
puted the polarization for solar wind ion velocities from 200- 
2200 km s"', using the cross sections published by Harel et al. 



al. |2005 1 assumes simple trajecto- 
and it is appropriate for describing 



( 1998 1. The upper limit of the considered velocity range cor- 
responds to ion velocities fou nd in the f astest Coronal Mass 
Ejections (CMEs) (Gopals wamy et al.|2009j , and the selected 
ion velocity range inclu des both the slo w (400 km s~') and 
fast (750 km s"') SW ( jSmifli et al.|2003| ). 

Two different models of the heliospheric plasma wer e used 
in our calculatio n. The classical "h ot model" (CHM) ( |LaIle- 
,ment 2004; Lallement et ; ' 
ries m a hot environment ; 
density distributions of hydrogen and helium averaged over 
short-term solar activity, including solar minima and maxima. 
This model has been employed in calculations of the helio- 
spheric CX X-ray background averaged on the solar activity 
cycle (Pepino et al. 2004). The CHM model does not include 
the termination shock and the distribution of charged SW par- 
ticles was assumed to be isotropic with respect to the Sun and 

to follow Nxe^(r) = N^-^l {Rq/k) . Since the heliopause is not 
defined within the CHM model, the integration limit was set 

to 300 AU. 

The second model ( |Muller et al.|2006[ ) is based on a MHD 
multifluid calculation with four interpenetrating fluids, where 
one 'fluid' represents protons of the interstellar and SW 
plasma, while three remaining components are used to model 
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Figure 1. Angular dependence of the polarization of CX X-rays emitted in a single de-excitation from the Ap (solid lines), 5p (dashed lines) and dp (dotted 
lines) excited states. Polarization for the angle of observation a is expressed as the radius-vector P(a) (%). The axis for o = 7r/2 (LOS in the upwind direction) 
is indicated, and the angle a increases counterclockwise. Blue cones indicate regions within ±10° of the Sun (right side) or directly opposite (left). Left panel 
shows polarization of the X-rays emitted by 0**'+\ calculated for the slow SW using the CHM (red lines) and MHD (black lines) models. Right panel shows 
polarization for the CX emissions from C*'^^' for the fast (green lines) and slow (black lines) SW calculated within the MHD model. 



the neutrals (see |Zank et al.| ( |l996 i for details). We used 
heliospheric hydrogen density distributions from the MHD 
model, and density distributions of C^^ and O^""" ions were ob- 
tained by scaling the proton density from the model accord- 
ing to i on abundances of the fast and slow SW ( [Schwadron &| 
Cravens 2000; Koutroumpa et al. 2006). The intensities were 



integrated up to the heliopause and plasma density variation 
across the termination shock was included in the calculation. 

For both models, we considered solar conditions corre- 
sponding to a typical solar minimum. All external pertur- 
bations that could mix populations of ionic excited states, 
such as heliospheric and interstellar magnetic fields, were ne- 
glected. The intensities were calculated in the entire ecliptic 
plane. 

4. RESULTS 

We calculated polarization maps of CX X-ray emissions in 
the ecUptic plane for C**^"^' and O**^"^* ions (Figure [T]i. The 
angular dependence of the polarization is a result of interplay 
of several factors. The parallel component of the intensity, /|| , 
depends on the angle of observation according to the factor 
( 1 - 3 sin^ 61) (in Eq. (pjl), while the total intensity is dependent 
on th e density of helios pheric plasma, which is rather com- 
plex (I Miiller et al.|2006 1. Namely, the ion density distribution 
is the greatest outside of the Solar corona and diminishes with 
the distance, roughly following the inverse square law inside 
of the termination shock. Neutral H is the dominant compo- 



nent of the interstellar flow. It penetrates into the heliosphere 
from (A,/3)h e = (252.3°, 8.5°) dire ction in the helio-ecliptic 
coordinates ( |Lallement et al.|2005| l, where H atoms are ion- 
ized as they encounter the solar wind and radiation. Conse- 
quently, the H gas density decreases as the ionization rate in- 
creases in the vicinity of Sun. The resulting H distribution is 
anisotropic and sensitive to the level of solar activity. Thus, 
the CX X-ray polarization map is affected by both the geo- 
metrical factor and the anisotropicity of the plasma density. 

We first computed the polarization maps of CX X-ray emis- 
sions from the Ap, 5p and 6p states of O**^"^' induced by the 
slow SW, where we used the CHB and MHD models to de- 
scribe the densities of H and ions for conditions typical of a 
solar minimum (Figure [T] left panel). Both models predict 
the highest polarization in the upwind direction of observa- 
tion that vanishes as the LOS angle approaches n. Ratios of 
polarization calculated for different initial states remain the 
same within the models, although the values calculated using 
the MHD model are about a quarter higher. The highest polar- 
ization is predicted for a = 68.3°, and a = 71°, for the MHD 
and CHM model, respectively. Quantitatively, polarization of 
the X-ray emissions from the 4p and 6p states in the MHD 
model is about P = 4%, and P < 3% for the initial 5/:' state. 
For comparison, the simple analytic model given in Eq. ([8]), 
predicts f = 13.4 % in the LOS direction a = 51.4° for erms- 
sions induced from the 4-p state by the slow SW, and scales 
similarly to the MHD model for the 5p and 6p states. 
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To better illustrate the dependence of the CX-induced X-ray 
polarization on the velocity, we show an analogous polariza- 
tion map for the slow and fast SW (Figure[r| right panel) con- 
structed for the MHD model. We predict lower polarization of 
the CX X-rays induced by the fast SW than by the slow SW. 
The only exception are CX emissions from C**^"'"'(5/?), which 
we find to have the polarization of about 7% for the optimal 
angle of observation, a w 70°. Although the fast SW is not 
present at low (below /3he = ±20°) heliographic latitudes, the 
calculated velocity-dependence of polarization may give ad- 
ditional insight in understanding the CX emissions induced at 
higher heliographic latitudes, as well as in similar astrophysi- 
cal environments. 

In Figure[2]we illustrate the polarization for a broader range 
of ion velocities for the three excited states of C**^"""' and O**^"^' 
in the upwind {a = 90°) and downwind {a = 270°) directions. 
The polarization for the optimal LOS, a = 68.3°, is about 17% 
higher than in the upwind direction. As an overall trend, the 
polarization increases as the ion velocity increases for ion ve- 
locities higher than 1000 km s"', except for the 6p state in 
C**^^"*"' where the polarization decreases for velocities higher 
than 1700 km s"'. This behavior is a reflection of dependence 
of the CX partial cross sections on velocity. While for high ve- 
locities we expect the populations of different magnetic sub- 
levels to asymptotically approach their statistical values, SW 
ion velocities are not high enough to clearly show that trend. 

5. CONCLUSION AND FUTURE DIRECTIONS 

In this work, we present the first calculations of polarization 
of the heliospheric X-rays induced by the CX collisions with 
fully stripped oxygen and carbon ions, O**"^ and C''"^. To ana- 
lyze its velocity and directional dependence, we calculated the 
polarization in the ecliptic plane as a function of the LOS for 
ion velocities from 200-2200 km s"', where the heliospheric 
plasm a was desc ribed usi ng two mo dels of different complex- 
ity (L aUement[200 4; M iiUer et al.|2000 ). Our calculations in- 
dicate that the CX heliospheric A-rays are mildly polarized 
and that the polarization depends on the SW ion velocity. 

While this study was restricted to C^"*" and O^^ ions in a 
simplified geometry, it nevertheless, illustrates a rather gen- 
eral property of the CX-induced radiation. In the considered 



interval of ion velocities colliding ion and atom form an inter- 
mediate quasi-molecule whose projection of electronic angu- 
lar momentum is quantized along the direction of the quasi- 
molecular axis. In collisions, quasi-molecular axes of differ- 
ent ion-atom pairs rotate until, at the end of encounter, they 
are oriented along the ion velocity vector. Thus, CX colli- 
sions yield an ensemble of aligned excited ions that emit po- 
larized photons. The polarization reflects a directional orien- 
tation of local collisional velocities in the regions of the SW 
plasma that are mostly responsible for production of the X- 
ray flux. Analogously, we may expect a relatively high level 
of polarization of the CX X-ray emissions from cometary at- 
mospheres or in Jovian polar X-ray auroras, where the ion 
velocity co-orientation in regions of the X-ray-producing CX 
collisions may be very high. Fast CMEs, propagating in the 
interplanetary gas, could also produce highly polarized CX 
X-ray emissions, particularly if the X-ray detector is located 
close to the CME's plasma stream. 

Polarization measurements could be used as a powerful 
tool to supplement the CX spectra and provide additional in- 
sight into underlying astrophysical processes. An example 
of such an environment, directly related to this work, is the 
Local Bubble. Several recent studies backed by new observa- 
tional data questioned the validity of the current Local Bub- 
ble picture, arguing that the contribution of heliospheric X- 
rays to the soft X-ray background may be higher than previ- 
ously thought (Robertson et al.'200r;'Koutroumpa et al."2007l 
Robertson et al. 2008 ; Henley & Shelton 2008 ; Welsh & SheL] 
ton|200 9 ). The polarization measurements could help identify 
the contribution from the heliospheric CX radiation in the dif- 
fuse X-ray background, especially if performed for different 
levels of solar activity. 

As a continuation of this work, a detailed map of the polar- 
ization of the SW induced X-rays in the Solar System should 
be constructed using results of an elaborate 3D MHD model 
of the heliospheric plasma and neutral gas. If computed for 
different solar conditions, such maps could be used as an in- 
dependent method for determining velocity and spatial dis- 
tributions of the solar wind plasma. Finally, the investiga- 
tion of CX-induced X-ray polarization may be extended to 
other interesting topics, including the He focusing cone, de- 
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pendence of polarization o n the solar activit y, or regions of 
the SW plasma turbulence ( Borovikov et al.|2 008 ). 
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